Electron properties of carbon nanotubes in a periodic potential 
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A periodic potential applied to a nanotube is shown to lock electrons into incompressible states that can form a 
devil's staircase. Electron interactions result in spectral gaps when the electron density (relative to a half-filled 
Carbon vr— band) is a rational number per potential period, in contrast to the single-particle case where only 
the integer-density gaps are allowed. When electrons are weakly bound to the potential, incompressible states 
arise due to Bragg diffraction in the Luttinger liquid. Charge gaps are enhanced due to quantum fluctuations, 
whereas neutral excitations are governed by an effective SU(4)~0(6) Gross-Neveu Lagrangian. In the opposite 
limit of the tightly bound electrons, effects of exchange are unimportant, and the system behaves as a single 
fermion mode that represents a Wigner crystal pinned by the external potential, with the gaps dominated by 
the Coulomb repulsion. The phase diagram is drawn using the effective spinless Dirac Hamiltonian derived in 
this limit. Incompressible states can be detected in the adiabatic transport setup realized by a slowly moving 
potential wave, with electron interactions providing the possibility of pumping of a fraction of an electron per 
cycle (equivalently, in pumping at a fraction of the base frequency). 

PACS numbers: 71.10.Pm, 85.35.Kt, 64.70.Rh 



I. INTRODUCTION 

Since their discovery,' carbon nanotubes (NTs) remain in 
focus of both basic and applied research. ^'^ Besides their 
important technological potential,"*'^ nanotubes are a testing 
ground for novel physical phenomena involving strong elec- 
tron interactions. In particular, they are believed to be perfect 
systems to study Tomonaga-Luttinger liquid effects.^" Ex- 
perimentally, effects of electron-electron interactions in nan- 
otubes have been observed in the Coulomb blockade peaks in 
transport, '^''^ in the power law temperature and bias depen- 
dence of the tunneling conductance, ''^'^ and in the power law 
dependence of the angle-integrated photoemission spectra.'^ 

Nanotubes are rich systems to study electron correlations 
for a number of reasons. Indeed, their one-dimensional (Id) 
character increases effects of interactions; spin and Brillouin 
zone degeneracy^ result in the presence of the four polariza- 
tions of Dirac fermions which can allow one to probe SU(4) 



spin excitations; 
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exceptional chemical and mechanical 



NT properties result in very low disorder; finally, diverse 
methods of nanotube synthesis, a variety of available nanotube 
chiralities and of the ways of coupling to the nanotube elec- 
tron system allow one to explore a wide region of parameter 
space. 

In the present work we propose to employ the coupling 
of an external periodic potential to the nanotube electronic 
system as a probe of both Tomonaga-Luttinger correlations, 
and of the Id Wigner crystallization, via commensurability ef- 
fects. We focus on the electron properties of single-wall NTs 
in a periodic potential whose period Acxt is much greater than 
the NT radius a, Acxt S> a. Such a potential can be realized 
using optical methods, by gating, or by an acoustic field. In 
all of these cases, the realistic period Acxt is of the order 0.1- 
1 /im. As shown below, effects of the Tomonaga-Luttinger 
coiTelations on the Bragg diffraction of electrons, realization 
of the SU(4) spin excitations, as well as pinning of the Wigner 
crystal can be demonstrated in such a setup depending on the 
applied potential and on the NT parameters. 



In particular, below we will identify incompressible elec- 
tron states, characterized by excitation gaps, that can arise 
when the average NT electron number density p (counted 
from half-filling) is commensurate with the period Aext of the 
external potential: 



TOtot 
Acxt 



TOtot = 4m . 



(1) 



In Eq. (1), m is the number of fermions of each of the four 
polarizations (herein called "flavors") per period. 

At what density values ni can the spectral gaps open? The 
single-particle treatment^^ maps the problem onto that of the 
Bloch electron, resulting in the spectrum of minibands sepa- 
rated by minigaps as a consequence of the Bragg diffraction 
on the external potential. Filling up an integer number m of 
minibands corresponds to adding m electrons of each flavor 
per "unit cell" (the period Acxt)- Thus, minigaps open up 
when the density (1) is integer, m — 0, ±1, ±2, ... In 
other words, the wave nature of the Bloch electron leads to 
commensurate states with integer density. 

Our main result here is that electron interactions dramat- 
ically change the spectrum, adding incompressible states at 
rational densities m = p/q, in which commensurability 
is induced by interactions.^' In a fractional— to state, the 
NT electron system is locked by the external potential into 
a (/Acxt— periodic commensurate configuration (such as the 
ones schematically represented in Fig. 1). Naturally, the states 
with the lower denominator q are more pronounced. Realis- 
tically, due to finite NT length and temperature, only a few 
states with small enough q can be detected. However, these 
fractional-TO states are important as the corresponding mini- 
gaps are interaction-induced (and vanish in the noninteracting 
limit). Measurement of such minigaps can provide a direct 
probe of interactions between electrons. 

From a practical standpoint, measurement of the gaps can 
identify the incompressible states characterized by quantum 
coherence on a macroscopic scale. A challenging experimen- 
tal proposal^" is to realize the Thouless pump^' by taking ad- 
vantage of the semi-metallic NT dispersion. In such a setup. 
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FIG. 1: (Color online) Incompressible electron states (schematic), 
(a) The (1/2,0) state with mtot = 1/2. (b) The (3,2) state with 
rritot ~ 1. Large potential amplitude breaks the "Dirac vacuum" 
placing holes (open circles) into the potential maxima, (c) An exam- 
ple of the phase soliton excitation over the (1,0) state. The excitation 
gap is given by the phase soliton energy, (d) The m = 1 incompress- 
ible state allowed by the Bloch theory of noninteracting nanotube 
electrons. Four kinds of labels mark the NT fermions of the four 
flavors, (e) Semiclassical picture of the simplest interaction-induced 
m = 1/2 state. Fermions of the same flavor avoid each other due to 
the Pauli principle. 

a quantized current is predicted to arise whenever the chemi- 
cal potential is inside the minigap created by the adiabatically 
slowly moving potential wave. As our approach suggests, 
due to electron interactions, commensurability (1) will result 
in additional current plateaus^' corresponding to pumping on 
average of a fraction of unit charge per cycle. Equivalently, 
electron interactions could allow one to realize a novel effect 
of adiabatic pumping of charge at the fraction of the base fre- 
quency of the potential modulation. 

In connection with the adiabatic current quantization, we 
note a parallel between our system and the quantum Hall ef- 
fect. In both cases, commensurability with external electro- 
magnetic field yields incompressible states, spectral gaps, and 
quantization in transport. Moreover, it is the electron interac- 
tions that result in the incompressible states at fractional fill- 
ing factors (commensuration due to interactions) in addition 
to the integer plateaus (commensuration due to wave nature of 
electrons), both causing non-dissipative transport.^^ This sim- 
ilarity extends onto the disorder stabilizing plateaus at lower 
denominators. 

From a theoretical viewpoint, the excitation spec- 
trum is linked to the general theory of commensurate- 
incommensurate transitions. ^^^^ In this approach, an excita- 
tion over the commensurate state (an incommensuration) is 
represented by a phase soliton [illustrated in Fig. 1(c)], whose 
energy gives the corresponding spectral gap. In the present 
work the phase soliton method is generalized onto the case of 



strongly interacting massive Dirac fermions of multiple polar- 
izations. 

Technically, the problem at hand requires non-perturbative 
treatment of electron interactions at low density (pAcxt = 
™tot ^ 1), i-e- at the bottom of the band, taking into account 
the effects of the curvature of the electronic dispersion. This 
is hard to achieve in the conventional bosonization scheme 
which rests on the linearized dispersion and describes hydro- 
dynamic density modes extended over the whole system. In- 
deed, the curvature of the dispersion is a "dangerously irrele- 
vant perturbation" that causes ultraviolet divergencies. How- 
ever, it is the curvature of the dispersion that can provide a 
length scale to an otherwise scale-invariant Gaussian bosonic 
action. Not surprisingly, phenomena in which the physics 
on the scale of Fermi wavelength is important, such as the 
Coulomb drag between quantum wires, often rely on the 
curvature and, hence, are difficult to treat. 

Nanotubes provide a crucial theoretical advantage: Even 
massive interacting Dirac fermions can be bosonized, by 
virtue of the massive Thirring — sine-Gordon duality, 
with the curvature of the dispersion controlled by the Dirac 
mass (NT gap Aq at half-filling). In the bosonic language, the 
gap controls the strength of the nonlinear term in the bosonic 
action. Moreover, the NT Hamiltonian is SU(4) symmetric 
in the forward scattering approximation. This symmetry sim- 
plifies the treatment since in the problem of the four coupled 
modes, only two different length scales appear: The charge 
scale, Ich (screening length for the Coulomb interaction), and 
the flavor scale, 1^ ^ Zch, with a meaning of a size of the elec- 
tronic wave function represented by a composite sine-Gordon 
soliton of the charge and flavor modes. '° 

External periodic potential locks electrons into incompress- 
ible states at commensurate densities. As a semiclassical illus- 
tration of this phenomenon, in Fig. l(d,e) the centers of soli- 
tons that represent electrons of the different flavors are marked 
by the four different symbols. In accord with the above, such 
a locking is only possible when the NT gap Aq ^ 0, since 
the dispersion curvature gives rise to the finite soliton sizes 

Zch, ^fl c>c Aq ''^^ with 7 > 0. For strictly linear dispersion (as 
in the armchair tube), both of these scales are infinite, and the 
minigaps do not open, as it is impossible to pin a system that 
has an infinite length scale. 

In general, both Bragg diffraction and the Coulomb interac- 
tion contribute to the excitation gaps. While the diffraction is a 
dominant factor for the incompressible states when the system 
is controlled by the Luttinger-liquid fixed point, the Coulomb 
interaction plays a major role when the electron system is in 
the Wigner-crystal regime. The external potential, by bring- 
ing the additional length scale Acxt to the problem, naturally 
distinguishes between the two sides of the Luttinger-liquid - 
Wigner-crystal crossover. Technically, this distinction occurs 
on the level of a saddle point of the nonlinear bosonic action. 

In a recent Letter,^' we showed that, in the semiclassical 
(strong coupling) limit, when electrons are nearly pointlike 
(Zfi ^ Acxt), the long range Coulomb interaction leads to 
the devil's staircase of commensurate states representing the 
pinned Wigner crystal. Gaps open at rational densities mtot, 
and Coulomb interaction sets the gap energy scale. An exam- 
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pie of such a state with mtot = 1/2 is shown schematically in 
Fig. 1(a). In this work we will study the phase diagram in this 
limit in detail. Remarkably, the Dirac character of nanotube 
electrons brings about a set of incompressible states in which 
the "Dirac vacuum" is broken when the potential amplitude 
exceeds the gap at half-filling. In this case, physically dif- 
ferent incompressible states can correspond to the same total 
density (1) with 

4™ = mtot = Ue -Uh- (2) 

To further characterize these states we specify the pair of num- 
bers {rig, rih) of electrons and holes in the potential minima 
and maxima correspondingly. An example of the (3, 2) state 
with mtot = 1 is schematically shown in Fig. 1(b). 

The opposite limit of the weakly coupled electrons (/a » 
Acxt) is connected to the gaps opening due to the Bragg 
diffraction. In this case electrons are delocalized over many 
periods, and gaps occupy a small part of the spectrum. Ex- 
change is important, and adds to the cost of an excitation in 
which an extra electron of a particular flavor is added to the 
commensurate configuration. Below we will specifically fo- 
cus on the lowest-denominator fractional state m = 1/2, for 
which we find the effective action that describes excitations 
over the commensurate state, and obtain the charge and the 
SU(4)-flavor excitation gaps. 

This paper is organized as follows. In Section II we intro- 
duce the model Hamiltonian. In Section III, as a warm-up, 
we consider the noninteracting problem. In Section IV the 
many-body Hamiltonian is bosonized. 

Section V is central, as we introduce the phase soliton 
method for the bosonized NT electrons. In Section VI we con- 
sider the Bragg diffraction limit, whereas in Section VII we 
focus on the pinned Wigner crystal. In Section VIII we out- 
line the phase diagram in the semiclassical (Wigner crystal) 
linnit. In Section IX we discuss experimental ways to detect 
incompressible states. 



II. THE MODEL 

Near half-filling, the nanotube electron system in the for- 
ward scattering approximation^'** is described in terms of the 
four (2spin X 2vaiiey) Dirac fermion flavors, with the following 
second-quantized Hamiltonian, 

H = 'Hn + Ti-hs + Ti-cxt ■ (3) 



The first term Ho is the massless Dirac Hamiltonian that in- 
cludes Coulomb repulsion: 




is a two component Weyl spinor of the flavor a, x is the 
coordinate along the tube, cri 3 are the Pauli matrices, and 
w « 8 • 10 cm/s is the NT Fermi velocity. In the total number 
density pk = Jdx e~''^^p(a;), with 

4 

p{x) ^^^i {x)^oc (x) , (6) 

a=l 

we neglected the strongly oscillating components ~ 
^±2tKBzx^ with Kbz ^ l/flcc defining the position of 
the Dirac points in the Brillouin zone of graphene [Ucc = 
0.143 nm being the length of the Carbon bond]. For a nan- 
otube of radius a placed on a substrate with the dielectric con- 
stant £, the Id Coulomb interaction 

V{k) = ^Vo{k), (7) 

where 

yo(fc)^e2ln[l + (fca)-2] . (8) 

The essence of the forward scattering approximation em- 
ployed above is that even in the presence of electron interac- 
tions, the two Dirac points at the corners of the Brillouin zone 
remain decoupled. Since the long range potential (7) does not 
discriminate between the Carbon sublattices, scattering ampli- 
tudes for the fermions of same and opposite chiralities at each 
Dirac point are equal. Hence in the Hamiltonian (4) we in- 
clude only the part of electron interaction that involves the 
smooth part (6) of the electron density. Furthermore, in the 
Hamiltonian (4) the backscattering and Umklapp processes 
between the Dirac points are discarded. These processes are 
parametrically reduced by acc/a ^ 0.1, while the Umklapp 
amplitude is also numerically small.' 

The second term in the Hamiltonian (3) describes the 
backscattering between the left and right moving fermions 
within each Dirac point: 

Hbs = ^0 y Ea=iV'^criV'aC^a; . (9) 

This backscattering is different from the usual "V{2kF)" term 
in quantum wires, as the interaction-induced gaps are unde- 
tectably small in metallic NTs. Instead, we rely on the cur- 
vature of the dispersion arising from the bare gap Aq at half- 
filling. Remarkably, this gap can greatly vary depending on 
NT chirality or external fields.^ In particular, semiconducting 
nanotubes have a large gap 

A(semicond) nioTr/ /in\ 

3^ eV/o[„m] (10) 

that scales inversely with the radius a. Metallic NTs can be 
of the two kinds. There are truly metallic, or the so-called 
"armchair" NTs which have a zero gap at half-filUng, Aq = 0. 
However, a gap Aq D that is small compared to the Id 
bandwidth 

Tiv 

£> = _ = 0.53 eV/a[,„] (11) 
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can appear due to the curvature of the 2d graphene sheet in the 
nominally metallic tube. This gap is inversely proportional to 
the square of the NT radius, and is numerically given by 

^(.emimet) ^ . j SOch | /af„^] (12) 

as a function of the NT chiral angle 8ch-^^'^^ Even smaller Ag 
can be induced in a strictly metallic "armchair" (Qch = 7r/6) 
tube by applying magnetic field parallel to the NT axis.^^'^^ 

Finally, interaction with the external periodic potential is 
represented by 

Wext = jdxpU{x). (13) 

In the present work, for the purpose of simpler algebra, we 
consider the potential of the form 

U{x) = A cos fcexta; , fcext = 27r/Aext • (14) 

Qualitatively, the results of this work will be valid for any pe- 
riodic potential realized by the means described in the Intro- 
duction. As the values of the period Aext are in the 0.1 — 1 fim 
range, the separation of scales fcext 1/a ^ l/occ rules out 
a possibility of coupling between the Dirac points via the po- 
tential (14), in full compatibility with the forward scattering 
approximation. 

In general, our main conclusion about the interaction- 
induced incompressible states at fractional densities that open 
up in addition to the integer-density gaps is valid for any value 
of the NT gap Aq. Quantitatively, our present analysis that 
will be based on the bosonization demands the separation of 
scales Aq <C D, practically requiring the semimetallic tubes. 

The Hamiltonian (3) written in the forward scattering ap- 
proximation is SU(4) invariant with respect to rotations in the 
space of the fermion flavors V-o - with both the Coulomb inter- 
action and the coupling to the smooth external potential (14) 
preserving this symmetry. 



in. NONINTERACTING ELECTRONS 

In the absence of interactions, the contributions of the four 
fermion flavors factorize. Here we consider the single electron 
spectrum [cf. Eqs. (4), (9), and (13)] 

{nD + U{x)}iP = ei}, (15) 

where Hd for fermions of each flavor is given by 

Hd = -ihvaadx + Aoai . (16) 

This spectrum has been briefly analyzed in Ref . 20. Below we 
will study it in detail emphasizing the distinction between the 
two opposite regimes of coupling, tight-binding vs. nearly- 
free electrons. It is convenient to perform a gauge transforma- 
tion 

^P'{x) = e-'^'^^^'"'=«*^ ^{x) , (17) 



where the dimensionless amplitude 

A = A/eo . (18) 

The kinetic energy scale 

eo = hkextv = /iw/Aext = 3.3 meV/Aext [nm] (19) 

is the Dirac level spacing in each potential minimum. 

After the gauge transformation (17) the Hamiltonian -|- 

U{x) becomes 

n' = -ihv d^az + AoCrie-2'-4^3 An^^.x ^20) 

Remarkably, due to the Dirac character of the problem, in the 
Hamiltonian (20) the relative importance of the potential en- 
ergy (the second term) is governed by the value of the Dirac 
gap Aq, rather than the potential amplitude A, whereas the ki- 
netic energy (the first term) is O{eo). [Note that in Eq. (20) we 
discarded the Schwinger anomaly term oc JdxU'^ (x) since 
its effect^^ of adding a constant energy shift is not important 
here.] 

The coupling of the NT electron to the external potential 
can be either weak or strong, depending on the gap value Aq. 
Consider the eigenvalue problem 

W'V' = e^' (21) 

that is periodic in Aext • Its solutions are the spinor Bloch states 
ip'pix) = Up(a;)e*''^ with a quasimomentum Tip taking val- 
ues in the effective Brillouin zone defined by the potential pe- 
riod, — fcext/2 < p < fcext/2. Below we distinguish between 
the weak coupling limit (nearly free electrons), and the strong 
coupling (tight-binding) limit. 
In the limit of nearly free electrons, 

i>p{x) = Up(a;)e'^^ , Up{x) ~ const, (22) 

and the electron wave functions are close to plane waves. Note 
that in the massless case of Aq = 0, Up = const, and the 
external potential has no effect since it is gauged away. It is 
a manifestation of the fact that the scalar external potential 
does not mix massless Dirac branches whose wave functions 
have a spinor structure. This is also consistent with the above 
discussion that the potential energy scale is determined by the 
backscattering Aq. 

When the backscattering Aq ^ 0, Bragg diffraction on the 
potential results in mixing between the left and right moving 
states of the Dirac spectrum at the values 

TTl 

Pm = ±-^K^t, m = ±l, ±2, ... (23) 

of the electron quasimomentum, opening minigaps in the sin- 
gle particle spectrum at energies = TOeo/2. Perturbation 
theory in Aq yields the minigap values^" 2Am\ 

A(^)(A)- Ao|J™(2A)|. (24) 

The superscript (0) in Eq. (24) relates to the noninteracting 
case. Minigaps (24) oscillate as a function of the potential am- 
plitude A and vanish for particular values of A corresponding 
to zeroes of the Bessel functions Jm- 
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FIG. 2: (Color online) Single electron spectrum in the potential (14) 
for the case Ao — eo- 



The opposite regime is the strong coupling, or the tight- 
binding limit, in which case the external potential localizes 
the semi-classical Dirac electrons (holes) in its minima (max- 
ima), tunneling between adjacent potential wells is exponen- 
tially suppressed, and minigaps become larger than the sub- 
bands. In this regime the wave function Up{x) becomes lo- 
calized around the potential minimum, with its characteristic 
size 



^wf ^ Ac 



(25) 



The spectrum can be obtained by solving the problem (21) nu- 
merically via the transfer matrix. For an illustration, in Fig. 2 
we plot the spectrum for the borderline case Ao = eo in which 
the potential and kinetic energy in the Hamiltonian (20) are of 
the same order. Similar to the nearly free electron limit, the 
energy gaps oscillate as a function of A and vanish at certain 
amplitude values. The spectrum also has a characteristic Dirac 
e — > — e symmetry. Fig. 2 suggests that even for moderately 
strong backscattering Aq = eo minigaps can occupy most of 
the spectrum when A eg. 

In Appendix A we show that the limit (25) of exponentially 
suppressed interwell tunneling holds whenever 



Ao 



A 

< — < 
eo 



Ao 

eo 



Ao > eo , 



(26) 



where the exponent c = 3 for tunneling from an energy level 
close to the bottom of the potential (or close to the top for 
holes), and c = 2 for tunneling from a level far from the po- 
tential minimum or maximum. 

Increasing the potential amplitude A beyond Ag/efj^^ ac- 
tually enhances the tunneling amplitude between the particle 
and hole continua. This is reflected in Fig. 2 by wider mini- 
bands for ^ > eo. Indeed, an increase in A makes the tun- 
neling barrier through the Dirac gap 2 Ao shorter, and reduces 
the action under the barrier. The other way of understading 



the apparent delocalization of the wave functions for large A 
is to notice that large potential amplitude results in strong os- 
cillations of the potential term in Eq. (20). These oscillations 
effectively average this term to zero in the limit A — + oo, in 
which case the minigaps become small as A~^^^. The latter 
follows from the large-A asymptotic behavior of the Bessel 
functions for the weak coupling minigaps (24) that become 
valid in this limit. 



IV. BOSONIZATION 

The Hamiltonian (3) is bosonized^'^''"''^'^" by virtue of the 
massive Thirring — quantum sine-Gordon duality, by rep- 
resenting the fermionic operators (5) as nonlocal combina- 
tions of bose fields — (27ra)~^/^ e'®°. The conjugate 
momenta 



h 



(27) 



obey the canonical relations [nQ,(x), G/3(y)] 
—ihSa/i 6{x — y). The result is the Lagrangian 



C 



hv X ^ 

a = l 



1 2 A 

- {dtQc^f - {d.&cf - cos 26, 

V nva 





- -Jdx'p{x)V{x-x')p{x')-p{U{x)+ti}, (28) 

where the density p{x) = Pa = \ d^Qa, and p. 

is the chemical potential (p ~ at half-filling). The Gaussian 
part of the action is diagonalized by the combination of the 
unitary transformation '"-^^ 



(29) 



with the density (6) being a gradient of the charge mode 0", 

pix) = ^d,d\ (30) 
and of the gauge transformation [cf. Eq. (17)] 

ea{x)^e^{x)-^fK-\u-p)dx' (31) 

that leaves 6'" intact, and shifts 

0'^O'~^rK-HU-p)dx'. (32) 
Here the charge stiffness 



Trnv 



(33) 



Below we drop the (iiTelevant) logarithmic dependence of 
the stiffness Kg on the momentum, assuming a constant 
value K = K^^i/i^^^, where ~ h is the charge soliton 
size (screening length for the Coulomb interaction). Using 
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e^/hv ~ 2.7 we estimate ~ 40 for Ich ^ 1 /im for the 
stand-alone tube; A' ^ 10 for the tube placed on a substrate 
with a dielectric constant e — 10. 

The nonlinear part of the action (28) is transformed using 
the identity ^Li cos 2&a = 4:J^{0°, 6*"), where 

3 3 

((9°, 9") = COS 61° • Yl cos 6'" + sin 6° ■ JJ sin 9", (34) 



a=l 
, 10,20 



a=l 



yielding the Lagrangian' 



(35) 



hv r 1 



27r Iv 

hv 

+ 



27r I w 

o=l 



(36) 



-Cbs = -^jr(0O + 2/ifcexta; - 2isin fcexta;, e"). (37) 

TT 

The Lagrangian (35) describes one stiff charge mode 6^ and 
three soft flavor modes 9", nonlinearly coupled by 

4Ao ^ £Ao 
hva a? D ' 



50 = 



(38) 



where D is the bandwidth (11), and the dimensionless quanti 
ties 

r A 



Keo ' 



Keo 



(39) 



are introduced in a way similar to that of Eq. (18), with the 
energy scale eo defined in Eq. (19). The difference between 
Eqs. (39) and (18) is in the screening (by a factor of 1/K) 
of external fields U (x) and /x by the interacting NT system. 
When go = 0, the electron number density (1) parametrized 
by m varies continuously with jl as 



m = 2jl. 



(40) 



In this case the system is a compressible scale-invariant 
Tomonaga-Luttinger liquid of the four flavors regardless of the 
electron interaction strength.^^ The nonUnearterm (37) breaks 
charge-flavor separation by binding O"" with 9^. It can become 
relevant at commensurate densities, yielding incompressible 
states. 

Note that the Lagrangians (28) and (35) are SU(4) invari- 
ant, as is the original Hamiltonian H. This invariance is not 
explicit in the adapted notation. It will manifest itself below 
on the level of renormahzation. 



V. PHASE SOLITON METHOD 
A. Electron as a composite soliton 

Evaluation of excitation gaps requires the knowledge of how 
an added electron is represented in the language of the cou- 
pled bosonic fields (29). Here we first discuss this issue in a 



simpler case of a stand-alone tube, when no potential is ap- 
plied: U{x) = 0. We follow the framework of Levitov and 
Tsvelik,'" emphasizing the details that will be important in 
the rest of the work. 

The idea is that when interactions are strong {K ^ 1), the 
action (35) is dominated by the charge sector, which justi- 
fies optimization of (35) in a saddle point fashion, with K^^"^ 
playing the role of 7i~^. In particular, the soft neutral modes 
0" adjust themselves to provide an effective potential for the 
charged mode 9^ in such a way that the cost of adding charge 
is minimal. For the most part one can regard the charge mode 
as classical, since its quantum fluctuations are suppressed by 
ii'-V2 (and can be further included as corrections to the scal- 
ing laws as described in Sec. V E below), whereas the neutral 
fluctuations strongly renormalize the coupling (38), with its 
renormalized value g given by Eq. (45) below. 

As a result,'" an electron added to a half-filled nanotube 
is a composite soliton characterized by the two length scales: 
the flavor l^, and the charge ^ch ~ -f^^^^^fl, so that Ich '> la- 
The neutral modes add a particular flavor to the electron, 
and, naturally, flavor is bound to charge by the non-linear term 
(37). Adding an electron of a given polarization a amounts to 
changing 8a by tt. According to the transformation (29), in 
the charge-flavor basis this results in a configuration in which 
all the fields 9''\ 6'" change by ±7r/2. However, given the 
separation of scales Ich ^ ^fl, solitons of the flavor modes 9°' 
appear as sharp steps on the scale of Ich that "switch" right in 
the middle of the charged soliton. The Coulomb interaction 
between electrons corresponds to overlap of charge sohton 
tails within the screening length Ich, whereas the fermionic ex- 
change corresponds to overlap of the (shorter) flavor solitons. 
In the low-density regime pla <C 1, overlap of the charge- 
soUton tails maintans quasi-long-range order of a Id Wigner 
crystal, whereas the effects of exchange are exponentiaUy sup- 
pressed. 

Technically, the effect of the neutral modes is to provide 
Fourier harmonics quadruphng in the nonlinear term (37), 
jr(6'°,r) T{9^) ~ cos46'0. The optimization leading to 
the effective potential T{9^ ) is shown in Fig. 3. On the classi- 
cal level, minimizing T{9^ ,9"-) over the values 9'^ = no7r/2 
with integer Ua yields the potential 

Tci{e^) = min 6*°) 

{e»=„„7r/2} 

(41) 

for the charge mode (top panel of Fig. 3). The switching be- 
tween the branches of ± cos 9^ and ± sin 9^ occurs via the 
neutral fields changing by ±7r/2 right in the middle of the 
slow charge soliton, when 9^ has changed by 7r/4. One ob- 
tains the classical soliton of the charge mode of the form'" 



min {±cos^°, ±sin 



2 cos ^ tanh(M — x/l^^^) , a; < 



where 



;(cl) 
'ch 



f - 2 cos-1 tanh(u + x/l^^ ) ' ^ > 
= {KjgYl'^ and tanhu = cos7r/8. 



(42) 



Quantum fluctuations of the neutral sector renormalize the 
bare potential ~ fifo-^cK^*^) by providing it with the scaUng 
dimension 



7 = 8/5. 



(43) 
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FIG. 3: (Color online) Classical Tdid'^) (top) and renormalized 
T{6^) (bottom) effective potentials (shown in bold) for the charge 
mode in the case of A" 3> 1, Eqs. (41) and (44) 



As a result, the effective potential for the charge sector 

-^bs — (^) Kihargc , Khargc = 5-^(^°) : 

where the renormalized nonlinear coupling 

1 Mo^^ 



(44) 



(45) 



and the cutoff D defined in Eq. (11). The potential T{6^) is 
plotted in the lower panel of Fig. 3. The result (44) is justi- 
fied in the adiabatic approximation /ch S> la when the poten- 
tial (41) is slow on the scale a < ^ < on which quantum 
fluctuations of the fields 0° accumulate (for renormalization 
procedure see Sec. V E below). 

Although the difference between the form of the potential 
(41) [considered in Ref. 10], and the renormalized one, T{6^) 
[Eq. (44)], is not crucial qualitatively,^^ it can be important for 
quantitative considerations. Namely, it affects both the soliton 
energy, its size, Zfl, and the effective screening length, Zch- Un- 
fortunately, the analytic form for the charge mode soliton with 
the potential term (44) is unavailable. One way to find this 
soliton is to approximate the potential (44) by a piecewise- 
quadratic function [this procedure works well'° for the soliton 
in the nonrenormalized potential (41)]. This yields 



9°(x) 



quad 



[x] 



4^ ' (46) 



x>0, 



with a charge soliton size 



Zeh = {Khgf^ , (47) 
which is a factor w 1.26 shorter compared to that in Eq. (42). 



Physically, strong Coulomb repulsion in the presence of the 
curvature Aq of electronic dispersion qualitatively changes 
the conventional "spin-flavor separated" Luttinger-liquid be- 
havior governed by the Gaussian action (36). At the Luttinger 
fixed point, all the modes enter on equal footing, and the dis- 
persion curvature is irrelevant at high density. At low density, 
the saddle point of the total action (35) describes the crossover 
to the Wigner crystal, in which flavor is bound to charge. The 
neutral and the charge sectors play very different roles in the 
Wigner crystal regime: the former creates the effective poten- 
tial for the latter This effective potential (44) has a period tt/2 
that is four times smaller than that of the potential (37) with 
fixed 9°-; the lowest harmonic in the Fourier expansion 



T{9°) = const - V /„ cos(4n6l°) 



(48) 



is cos 40°. The purpose of this period reduction is to lower 
the Coulomb energy by splitting the f?" — > 6*" + 27r excitation 
that carries charge 4e [according to Eq. (30)] and is a flavor 
singlet, into four subsequent excitations each carrying a single 
fermion of a unit charge and of a particular flavor. 

Practically, to obtain a smooth approximate form of the 
charge soliton, and to estimate its energy that dominates the 
cost of adding an electron, here we keep the most relevant 
harmonic in the Fourier series (48), 



jc-(6i°) w -/i cos 46*", /iw 0.178. 



(49) 



The approximation (49) overestimates the classical charge 
soliton energy by 10%: 



jTiappiox. 

£^chsoi. =ci — VKg, 



Cl 



0.844 . 



(50) 



compared with the energy in the full potential (44): 

^ch.sol. ~ 




(51) 



VS / d9° \/\ ~ cos-r d" w 0.752 . 



The charge soliton in the potential (49) takes the form 

^'approx(2;) = ^ cos'^ tanh(-2://ch' ) (52) 

with the rescaled size l^h = jiK/fig)^^^ w 'chVl-69. Note 
that, due to strong renormalization by the neutral fluctuations, 
the potential (44) is steeper than its classical counterpart (41). 
This leads to shorter charge soliton length than that of a clas- 
sical solution (42), as both of the approximations, Eqs. (46) 
and (52), indicate. 

Remarkably, the difference between the Luttinger liquid 
[described by the Gaussian Lagrangian (36)], and the Wigner 
crystal will become crucial after adding the external potential 
U (x). Depending on the relation between its period, Aoxt, and 
the "size of an electron", Ifi, either the Gaussian saddle point 
{go — > 0) of the action (35), or the Wigner-crystal saddle point 
described above will play out. The latter one applies in the 
regime Acxt ^ Ifi considered below in Section VII. There we 
show that the system effectively behaves as that of a single 
mode. 
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B. Incompressible states of the classical bosonic fields 

The above discussion of the transition to the Wigner crystal 
suggests that, already on the classical level, the nonlinear ac- 
tion (28) or (35) captures important physics even on the short 
length scale (smaller than the inter-particle distance!). This 
intuition can be readily extended onto the U ^ case. By 
adding the external potential (14), below we use the classical 
Hamiltonian 



K - 1 



■jgoJ2cOs2ea + ^da:Qa ' {u {x) ~ ^1 



(53) 



that follows from the Lagrangian (28), to numerically illus- 
trate how incompressible states appear For convenience here 
we work in the original basis 8q. 

Consider the simplest fractional density, m — 1/2, corre- 
sponding to the chemical potential jl — 1/A [cf. Eq. (40)]. 
Its classical ground state is an incompressible configuration in 
which the solitons of the fields occupy every other poten- 
tial minimum, as shown in Fig. 4. This Figure is a result of the 
numerical minimization of the Hamiltonian (53) with respect 
to the fields 8q. In agreement with the general theory (de- 
veloped for a single-mode system), for the fraction m = p/q 
the density should have the period qXcxt- In our case of the 
four modes, naturally, the period of the density ^dxQa in 
each mode is gAoxt, with q — 2 in Fig. 4. Note that, since 
the total density mtot = 2 is integer, the charge density pe- 
riod coincides with that of U{x). In the absence of U{x), 
all the solitons of 8 a would be equally separated from each 




FIG. 4: (Color online) Classical ground state of the Hamiltonian (53) 
for the density m = 1/2. Bold solid line marks the external poten- 
tial. Fine lines in the upper part denote the solitons of the fields 
Oa found by the numerical minimization of Eq. (53) for jl — 1/4, 
K = 10, go = 12fcext, A — Keo- The corresponding fine lines 
in the lower part denote the flavor densities Sp^ — Pa — rn/Xc^t 
counted from their average value m/Acxt. Bold dashed line is the to- 
tal density counted from its average value 4m / Aext . Note the period 
doubling for Qa{x). 



other due to the Coulomb repulsion. The finite— t/ configura- 
tion shown in Fig. 4 is the result of an interplay between the 
mutually repelling solitons and a confining periodic potential. 
We note that the fermionic exchange of the original problem 
(3) is manifest here akeady on the classical level, in the fact 
that solitons of the same flavor repel each other stronger, and, 
therefore, do not occupy the adjacent potential minima. 



C. Phase soli ton method for a single mode 

To systematically study excitation gaps over commensu- 
rate configurations, we employ and generalize the phase soli- 
ton method. Representing an excitation in a commensurate 
phase by a phase soliton was utilized in the past in various 
contexts. The model of locking a system into a commen- 
surate state was first suggested in the work of Frenkel and 
Kontorova.^"^ Later, it was re-discovered and solved by Frank 
and van der Merwe^^ in the context of atoms ordering on crys- 
tal surfaces, and by Dzyaloshinskii^^ describing a transition to 
the state with a helical magnetic structure. The general theory 
of commensurate-incommensurate phase transitions has been 
finalized by Pokrovsky and Talapov.^^ The result of these in- 
vestigations is a finite or infinite sequence of gaps ("devil's 
staircase") corresponding to commensurate states. 

The essense of the phase soliton method can be demon- 
strated on the example of the single-mode system. Consider 
the model Hamiltonian of the form 



■Hi 



- g cos{f39 + mkcy±x + 6 sin fcoxta;) (54) 



[cf. Eqs. (28), (35), or (53), albeit with only one mode 
present]. To make contact with the nanotube system, the pa- 
rameters in the Hamiltonian (54) have a similar meaning to 
those in the action (28) transformed by (31), with (3 = 2 and 
only 8i present. The Hamiltonian (54) is the continuous form 
of the Frenkel-Kontorova model^'* that describes locking of 
the compressible lattice of interacting particles (e.g. a chain 
of atoms) by the external potential of the period Aoxt (e.g. the 
substrate potential for the atomic chain), as is schematically 
shown in Fig. 1(a). The phase field represents the deviations 
from equilibrium positions along the chain, with changing 6 
by 27r//3 corresponding to a shift of the chain by its period. 
The gradient (38x9/2^ thus has a meaning of the excess par- 
ticle density near point x. The coupling between the particles 
and the potential is represented by g, whereas the lowest har- 
monic of the potential is cx 6 cos kcxtx- The crucial parameter 
is the period ratio m (number of particles per period Aoxt)- 

External potential can lock the chain into incompressible 
states when m is either an integer, or a simple fraction m = 
p/q. An excitation (an incommensuration) over such a com- 
mensurate state corresponds to adding an extra particle to the 
system. This changes the positions of other particles in the 
chain due to their mutual interaction. When the coupling g is 
small, such a change occurs over many periods Aoxt [as illus- 
trated in Fig. 1(c)]. 

In the language of the continuous model (54), the incom- 
mensuration is represented by the soliton-like change of 6 by 
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27r//3, with the excitation gap given by the soliton energy. 
The effective theory for the excitation follows from the Euler- 
Lagrange equation for (54), dl-Li/SO = 0, practically found 
by expanding in the powers of g, 

= + 0W + + 6»(") = 0(5"). (55) 

Here 9 is constant in the commensurate phase, acquiring a 
slow coordinate dependence 9{x) in an excitation {phase soli- 
ton). When m is integer, the zeroth order term = ^ in the 
expansion (55) is enough: Averaging over a period _Aext> the 
effective Hamiltonian is K (Oxd)"^ / 2 + g Jm{h) cos/36, and the 
excitation is described by the sine-Gordon soliton of a slow 
mode ^(a;) that interpolates by 27r//3 over a few Aext- 

For the fractional density m = p/q, averaging gcos{/39 + 
mkextx + b sin kc^tx) over qXcxt gives zero for constant 9. 
This indicates that the effective potential for 9 is of the higher 
order in the coupling g. Using the expansion (55), one obtains 
the lowest order Hamiltonian for the phase soliton 9{x) of the 
form K{dx9)'^ /2 + gq{b) cos /3q9, g^ ^ 9'. The correspond- 
ing soliton energy estimates the excitation gap. One also has 
to average over the thermal or quantum fluctuations to assess 
the relevance of the potential term ^ cos [3q9, in which case 
the (omitted) overall coefficient in front of the Hamiltonian 
(54) becomes important, as it controls the size of fluctuations. 



D. Phase soliton method for a nanotube 

Below we show that the case of multiple modes is qual- 
itatively different from a conventional single mode situation 
described above. This difference manifests itself in the pres- 
ence of the two regimes of coupling of the NT electron system 
to the external potential. Physically, the choice of the regime 
depends on whether the electron wavefunctions overlap with 
each other substantially, — that is, whether the physics is 
determined by the quantum interference with fermionic ex- 
change being important, or by the Coulomb repulsion between 
nearly point-like electrons, in which case the exchange effects 
are negligible. 

Technically, the difference between the regime of Bragg 
diffraction in a Luttinger liquid ("weak coupling"), or of pin- 
ning of the Wigner crystal ("strong coupling") is determined 
by the saddle point of the nonlinear action (35). In particular, 
if the neutral modes switch fast on the scale on which the ar- 
gument of the potential (37) changes appreciably, the Wigner- 
crystal saddle point (Sec. V A) is selected, and the system is 
described in terms of the single (charge) mode, with the prob- 
lem reduced to that described in Sec. V C above. If, on the 
other hand, the neutral modes are slow, the weak-coupling 
saddle point (expansion around go = 0) is the relevant one. 

To define the regimes of coupling we need to compare the 
size /fl of the phase soliton 9°'{x) of the flavor sector (size 
of the added electron) with the potential period Aext- For now, 
we assume that /a is known. The scale /a will be later obtained 
self-consistently as a result of integrating over the quantum 
fluctuations (performed for different cases in Sections V E, VI 
and VII). 



Weak coupling 

In the weak coupUng regime the flavor soliton tails, which 
correspond to fermions of the same flavor, overlap. Physi- 
cally, it means that the system "knows" that it is composed 
of the particles of different flavors since the role of exchange 
is important. An extreme example of this regime is a non- 
interacting system in which fermions of the same flavor ef- 
fectively repel due to the Pauli principle, whereas fermions of 
different flavors do not notice each other. 

Strong overlap between the soliton tails in a commensurate 
configuration requires large Ifl ^ Aext, or, equivalently, small 
coupling 5o in Eq. (35). The system is close to the Luttinger 
liquid, i.e. it is compressible for almost all densities apart from 
a few commensurate ones, for which small minigaps open up. 
This regime has been considered in Ref . 20 for integer m. The 
smallness of go warrants finding the effective Lagrangian 

,9'']^ Co [9" , 9-] + Cl"^^ , r ] (56) 

for the phase modes 9° and 9"' perturbatively in go (simi- 
larly to the single-mode model case considered in Sec. V C 
above), treating all the modes on equal footing. Namely, from 
the Lagrangian (35) with the chemical potential (40) corre- 
sponding to the commensurate density (1) parametrized by 
m = p/q, one finds the effective potential for 
the phase modes to the lowest order in 170 by solving the 
time-independent Euler-Lagrange equations = Oand 

5£'/69°' = 0, via the expansion 

r = + 6'"(^) + ... + + (57) 
= O(5o") , 9-^-^ = 0{g^) . 

The first term in (56) is the Gaussian Lagrangian (36) as a 
function of the slow phase modes, and the obtained potential 
energy £ JJJ* is of the order g^ (before integrating over quantum 
fluctuations). 

From the Lagrangian 9""] one finds the commensu- 

rate ground state in which the phase modes 9'^ and 9°' are con- 
stant. Similarly to the electron added to a stand-alone tube 
being represented by a composite soliton (Sec. V A), an ex- 
citation over such a ground state is a composite phase soliton 
of the slow phase modes 9°{x), 9"'{x). It can be found as a 
saddle point of the Lagrangian Cm , in a way similar to 
that described above in Sec. V A. 

The weak coupling limit is illustrated below by considering 
the simplest situation of integer density m, in which case the 
decomposition (57) is trivial since it contains just one term 
for each mode. Hence Cm [^°, 9°'] is precisely the Lagrangian 
(35) written as a function of 0°'" = 9^'"^ with the chemical 
potential jl given by Eq. (40). 

Technically, the weak coupling regime is characterized by 
the slow "switching" of the flavor modes 9°- on the scale on 
which the phase of the charge part of the potential energy 
C'^*[9°, 9"] in the Lagrangian (56) changes by 27r. For in- 
teger m, the potential energy is just the nonlinear term (37) 
written as a function of the slow modes 9°, 9°'. Thus we de- 
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mand 



2 A sin fccxt^n ^ Stt . 



(58) 



The first condition in Eq. (58) requires the flavor excitation to 
be extended on the scale of the separation 4p~^ = Aext/^i 
between same-flavor fermions. The meaning of the second 
condition will be made more clear below. The weak coupling 
limit (58) results in the following approximation for the en- 
ergy (37): 



(59) 



Here we discarded spatially oscillating terms (in other words, 
we averaged T over the period Aext), and denoted by the ar- 
row the soft mode "switching" that produced the optimized 
potential (41). 

The procedure (59) as written is allowed for the integer den- 
sity rn. In the case when the density m = p/qisa simple frac- 
tion one needs to utilize the phase sohton approach to find the 
effective potential of order and then perform analogous op- 
timization. The weak coupling limit is considered in Sec. VI A 
for integer m and in Sec. VI B for the simplest fractional den- 
sity m = 1/2. 



Strong coupling 

In the opposite, strong coupling limit, the Coulomb inter- 
action wins over the effects of the fermionic exchange. In this 
regime the system "forgets" its four-flavor nature, and effec- 
tively behaves as that of spinless Dirac fermions with the total 
density (1). Technically, the soft mode "switching" (denoted 
by the arrow below) occurs before averaging over the period 
of the potential, 

J" (e" + mfcexta; - 2isinfcexta;, 

f {W' -h mkextx — 2 A sin fcext^;^ 

= /i cos ^4^° -I- mtotfcexta; - 8 A sin fcexta;^ 

+/2 cos ^8^° -I- 2mtotfcexta; - 16i[ sin fcext^^) + ••• 

« /i cos (^46° + mtotfcexta; - SA sin fcext^;) • (60) 

Here the Fourier coefficients /„ are defined in Eq. (48), and 
in the last line we used the approximation (49). Note the de- 
pendence of the resulting potential energy on the total density 
mtot = 4m. The condition for the saddle point (60) is 



_ m/fl < Acxt , 
2^sinA:oxt^fi < Stt . 



(61) 



This regime is considered below in Section VII. 
Now let us clarify the meaning of the second condition in 

Eqs. (58) and (61). Increasing the potential amplitude A much 
beyond KHv/Ir effectively averages out the nonlinear term 



and thus draws the system into the weak coupling limit (58). 
This situation is similar to the noninteracting case considered 
in Section III above, where the approximation (24) becomes 
valid at large A. 

To summarize, above we have generalized the phase soliton 
method to the multiple-mode case. After the appropriate ef- 
fective action for the slow modes is chosen [either Eq. (56), 
or the one based on the potential energy (60)], excitations 
are described by the corresponding phase solitons. Due to 
the electron-hole symmetry of the original problem, adding 
an electron to a commensurate state (phase soUton) costs the 
same energy as removing an electron from the same state 
(an ti- soliton), while the sum of the two energies, 2Am, is the 
corresponding excitation gap. 



E. Effect of quantum fluctuations 

To illustrate integration over quantum fluctuations, below we 
renormalize the gap Aq for the stand-alone tube. (General- 
izations for the [/ ^ will be made in subsequent Sections.) 
Fluctuations of the neutral sector around the saddle point de- 
scribed in Sec. V A above are governed by the Lagrangian 
(35) with fixed value of 9°. Adiabaticity of the charged mode 
at K ^ 1 ensures that 9^ ~ const on the length scales 
a < I < la where quantum fluctuations of the flavor sector 
accumulate. With each neutral field contributing to the scal- 
ing dimension by 1/4, flavor fluctuations result in the scaling 
dimension 

70 = 3/4 (62) 
for the adiabatic charge mode potential (41), yielding 



Kharge(a) = 90^c\ 



Charge (0 — Charge (^) 

(63) 

where I is the renormalization-group (RG) scale, and a is the 
tube radius. Since 70 < 2, the potential Charge is relevant and 
grows. The flow (63) stops on the scale Z on which the 
potential energy becomes comparable to the kinetic energy of 
the flavor sector, since the perturbative RG yields the law (63) 
only while the renormalized coupling stays smaU.^' For the 
larger scales I > la the potential energy dominates and the 
problem becomes classical. 

The scale la has a twofold meaning. First, it is the correla- 
tion length for the flavor, estimated self-consistently from the 
balance of kinetic and potential terms 



Charge (^fl) — . 2 



(64) 



Beyond a; ~ Za the correlation functions of the fields de- 
cay exponentially, rather than in a power-law fashion. From 
Eqs. (63) and (64) we obtain the renormahzed potential (44) 
for the charge mode with the coupling (45) and the scaUng 
exponent (43), 



■ 70 



(65) 
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The second meaning of the scale Ifi is the renormaHzed ij'ze of 
the flavor soliton of the model (35) with fixed 0°. 

Since the charge sector is stiff, the excitation gap A is dom- 
inated by the charge soUton energy. The latter can be now 
estimated classically [cf. Eq. (50)] from the effective charge 
mode Lagrangian 



L 



charge 



charge 



where we used the saddle point approximation (44), (45) and 
(49) outlined in Sec. V A above. In the limit K ^ oo the 
effective Lagrangian (66) yields the gap"''^^ 



(67) 



At finite K the problem of finding exact scaling laws in the 
nonlinear theory (35) is far from trivial. Here we attempt to 
estimate the 0{r]) corrections, 



V 



(68) 



to the K = 00 scaling laws in Eq. (67). For that we note that 
in one-loop RG, one averages the nonlinear term £bs over the 
independent Gaussian fluctuations described by the Gaussian 
action Cq diagonal in charge and flavor. Due to the factoriza- 
tion (34), the result of such an averaging can be represented 
in the form of the product of the contributions of independent 
Gaussian theories, with the terms of the action (35) scaUng as 



K 



1 



-50 



hi 



-70 



-r,/4 



(69) 



Minimizing the composite-soliton action of the form (69) over 
both Ich and amounts to including the couphng between 
the flavor and charge fluctuations by the nonlinear term. The 
variational estimate yields 



a \Ao 



V 



10-2)7 



(70) 



for the correlation lengths, and produces the renormaUzed gap 
A ~ hvK/lch that scales as 



A~ ifs-CDE^Ao"- 



(71) 



The latter expression smoothly crosses over to the noninter- 
acting case -ftT = 1, A = Aq, and it also has a correct scahng 
(67) in the classical limit K ^ oo^^ 



VI. WEAK COUPLING LIMIT 

A. Integer density m = ±1, ±2, ... 

The saddle point (59) maps the problem of finding excita- 
tion gaps onto that with [7 = and 



90^ 9m- goJuiC^A) . 



(72) 



Based on the results of Sections V A and VE, one obtains the 
renormalized minigaps 2Am, where 

A„^/fV2-C^I^(AW(i))^ (73) 

[cf. Eq. (71)], and where the bare minigaps 

A<fi\A) = Ao \Jm{2A)\ (74) 

given by their noninteracting values (24) with the screened 
potential amphtude A defined in Eq. (39). Eq. (73) in the 
K ^ 00 Umit corresponds to minigaps found in Ref. 20. 
The qualitative features of the noninteracting minigaps (24) 
persist in the interacting case.^" Namely, as a function of the 
screened potential amplitude (39), the minigaps (73) oscillate, 
vanishing at particular values of A. However, minigaps (73) 
are strongly enhanced in magnitude compared to (24) due to 
electron interactions. Also the dependence of the minigaps on 
both the bare backscattering Aq and on the periodic potential 
amplitude A has a characteristic power law behavior which, in 
the limit of strong interactions K ^ 1, is given by a univer- 
sal power law 4/5 characteristic of the number of NT fermion 
polarizations at half-filling. 

What is the cost of adding electron's flavor on top of its 
charge? According to the mapping (72), the problem is 
formally equivalent to that with [/ = 0, so the results of 
Ref. 10 apply: On the energy scale below that of the frozen 
charge sector, the flavor sector is governed by the effective 
SU(4)~0(6) Gross-Neveu Lagrangian'"-''! 



ffGN^cos26'" cos 26''' 

a>b 



= hv Jdx{ixj^i,d^Xj - 9GN{XjXj){Xj'Xj')} ,05) 

where Xj are the Majorana fermions and the Gross-Neveu 
coupling in our case is gcN oc gm- The excitations of the 
model (75) are massive relativistic particles transforming ac- 
cording to different representations of the 0(6) group,"*^ with 
the mass scale physicaUy originating due to effects of ex- 
change. 

B. Fractional density m = i 

Using the phase soliton method outlined in Sec. V C above, 
we derive the effective phase mode Lagrangian 

£1/2 , r ] = Co [e" , + [^0 , r ] (76) 

given by the sum of the Gaussian part £0 [Eq. (36)] and the 
potential energy [see Appendix B for details] 



-4% = ^|(4-K)t;i/2(2^)^(20"°,2r) 

-KWl/2(2i) ^ cos Id" COS 26° 

a 

+KUy2{2A) ^ cos 2^" cos 26^ I . (77) 

a>b ) 
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Here the function T is defined in Eq. (34) and the other quan- 
tities are defined in Appendix B. 

Let us discuss the potential energy (77). Its first term has 
scaling dimension 3 and is irrelevant. The second term of 
the potential (77) is responsible for the charge excitation gap 
^ch,i/2- In the limit K ^ 1, this term has scaling dimension 

71/2 = 1 (78) 

and grows under the renormalization group flow. Finally, the 
third term is marginal. It describes the SU(4) flavor physics 
on the energy scale Ah ~ iir~^/^Ach,i/2- 

The charge gap is found via the saddle-point optimization 
of the second term of the potential (77) by the neutral sector 
in a way similar to that of Sec. V A above. Adding an extra 
electron to the system (76) corresponds to a composite phase 
soliton in which the 0" fields "switch" by 7r/2 in the middle 
of a slow charged phase soliton, at the point when cos ^ 
0. (It is possible to show^^ that although at this very point 
the effective sine-Gordon coupling cx cos 2^" for the neutral 
sector vanishes, the flavor soliton has a finite size and energy. 
This is the case since away from the center of the charged 
soliton, cos 2^° 7^ giving the finite flavor scale.) 

On the classical level, the neutral sector produces the charge 



potential ^j.^^^ [top panel of Fig. 5], 



^cos26'°cos2r ^ J'^^l^^\e°) = - |cos26'° 



(79) 



By integrating over the flavor fluctuations we obtain the flow 
of the form (63) with the scaling exponent (78) for the poten- 
tial energy gi/2^ci^^^ (^°), where 



51/2 = iKg^l^'i/2(2i)| . 



(80) 



From the self-consistency of the form (64) we find the flavor 
scale 



^fl,l/2 



£0 



K|t-i/2(2v4)| VAo. 
[a is the tube radius], the renormalized coupling 

1 



5i/2 - 

and the scaling exponent 



/2 



« gl/2 : 



7l/2 



(81) 



(82) 



(83) 



for the classical optimized potential (79). As a result, the ad- 
justment of the neutral sector yields the following effective 
potential for the charge mode: 



(1/2) 



charge 



JF(1/2) 



(1/2) 



7l/2 



■ cos4f?° + const. 



(84) 



Note that, due to the specific value of the scaling (78) [re- 
markably coinciding with the gap scaling for a single mode of 




FIG. 5: (Color online) Classical J-'^^'^\6^^) (top) and renormal- 
ized jt'^''^' (6'^) (bottom) effective potentials (shown in bold) for the 
charge mode at density m = 1/2 in the case of A' 1, Eqs. (79) and 
(84) correspondingly 



noninteracting fermions!], the renormalized potential JT^^/^) 
has only one Fourier harmonic, hence the coefficient 1/2 in 
Eq. (84) is not approximate [as one could anticipate from the 
analogous procedure leading to Eq. (49)], but exact. The func- 
tions .Fjfy^' and T^^l'^^ are shown in Fig. 5. They both have 
a period 7r/2 that corresponds to adding unit charge accord- 
ing to Eq. (30). As a result, the charge soliton for m = 1/2 
in the limit K \ has the exact form (52) with the size 



ch,l/2|if»l 



i(2i^/5j/2)'/'-^'/V.9i/2. 



For large but finite the soliton sizes /ch,i/25 ^VlsIi^ ^"d 



the renormalized minigap A 



ch,l/2 



ch,l/2 



found from the variational estimate of the form (69), 



K 



ch,l/2 



^fl 1 



/2 



.91/2 



'fl,l/2 



^ch.1/2 



can be 



(85) 



yielding \r] = i^^^/^, k ^ {K ~ 1)/K] 



Ach,i/2 D 



K\vif2i2A)\ 



Ao 



(86) 



The value of the gap (86) is strongly enhanced by the band- 
width D [even stronger than for the integer— to case of 
Sec. VIA] due to flavor fluctuations. We also note the to — 
1/2 incompressible state is explicitly interaction-induced. In- 
deed, although the charge excitation gap Ach,i/2 is derived in 
the strongly interacting limit K ^ 1, Eq. (86) gives a correct 
noninteracting limit Ach,i/2 = expected from the single 
particle Bloch theory. Formally, when K — 1, both flavor and 
charge gaps are zero since the last two terms of the potential 
(77) vanish (k 0). 

The weak coupling estimate (86) is valid when the flavor 
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soliton size (81) is large, 

^fl,l/2 > Aext , (87) 

according to the condition (58) above. Practically, due to 
the large bandwidth D, Eq. (87) requires rather small bare 
gap Aq. For typical parameter values, Aext ~ 0.1 /um and 
Ao = lOmeV, the soUton scale defined in Eq. (81) is small 
compared to the potential period Aext, 'fl,i/2 ^ 10 nm < Acxi, 
and the condition (87) does not hold. In this case the excita- 
tion gap is given by the strong coupling limit expression (103) 
(Section VII below). 

However, the result (86) is applicable for realistic param- 
eter values, whenever the flavor soliton size Zfl 1/2 becomes 
large: Either when one takes the tube that is almost metal- 
lic, Ao ~ ImeV, in which case Zfl,i/2 ^ 1/um > Aext. or 
close to certain potential amplitude values A* that correspond 
to zeroes of Vii2{2A) (Fig. 8). 

In the special case when the potential amplitude is such 
that the coupling vanishes, Wx/2(2^) = 0, the charged 
excitation is gapless [similar to the integer m case consid- 
ered above, where the gaps (73) vanish at the zeroes of the 
Bessel functions], but the flavor sector remains gapped. Its 
gap Afl can be estimated from the effective 0(6)~SU(4) 
Gross-Neveu Lagrangian of the form (75) that is given by 
the last term of the Lagrangian (77), where now the coupUng 
^GN oc ui/2{'2.A*) ^ (Fig. 8). We stress that the resulting 
Gross-Neveu coupling is a functional of the applied potential 
and thus, to an extent, can be controlled externally. Phys- 
ically, the origin of the flavor gap corresponds to the addi- 
tional cost of destroying the flavor order in a commensurate 
state by adding an extra electron of a particular flavor. Nat- 
urally, Afl ^ iiT^^/^Ach is controlled by the strength of the 
exchange interaction. 



VII. STRONG COUPLING LIMIT 
A. The Lagrangian for the charge sector 

Under the condition (61), the Wigner crystal saddle point 
(60) combined with integrating over the flavor fluctuations 
yields the following Lagrangian for the charged mode [cf. 
Eq. (66)]: 

-gfi cos(46'° + mtotkextx - 8iLsinA;exta^)| . (88) 

Here the density (1) is related to the chemical potential /i by 
Eq. (40), 8jl mtot. and the coupUng constant g is given by 
its renormaUzed value (45). 

The Lagrangian (88) describes the effective theory of a sin- 
gle mode in the presence of the external potential that is adi- 
abatic on the scale Zfl. This is a result of integrating out the 
three neutral modes, whose effect has been to produce the 
Wigner crystal saddle point [Sec. V A] and to renornalize the 



coupling g. As a result, the number of degrees of freedom in 
the theory (88) is significantly reduced. This simplification 
is justified under the following assumptions: (i) separation of 
charge and flavor energy scales Ach ^ Afl (for that, strong 
interactions -FT » 1 are crucial); (ii) dilute limit pla <C 1 that 
ensures exponentially small flavor-flavor correlations. These 
conditions imply that the flavor physics is decoupled from the 
charged one. In other words, the allowed flavor configura- 
tions in a train of composite solitons that represent electrons 
at low density enter with equal weights, which allows us to 
trace over them, effectively eliminating them from the action. 
Physically, such a situation corresponds to the limit when the 
temperature T ^ and exchange scale Afl ^ in such a 
waythat Afl/T^ 0.''3''*4 

We stress that in the limit (61), the renormaUzed coupUng 
g ~ l//fl^ in the Lagrangian (88) does not depend on the 
density and on the external potential. This is contrary to the 
weak coupUng case of Section VI where the renormaUzation 
of nonlinear couplings was sensitive to m and A. The reason 
for this sensitivity in the weak coupling limit is in the flavor 
soliton size la (that controls the scale over which the fluctua- 
tions are accumulated) extending over a few potential periods. 
Naturally, in that case the form of U{x) must influence the RG 
flow. 

In the present case, la <C Aext, and the RG flow produces 
the effective coupling g [Eq. (45)] on the scale a < I < la that 
appears microscopic for the external potential (14). In other 
words, the adiabatic (on the scales / < la) external potential 
has the effect of the chemical potential that can at best affect 
the local charge density, but not the RG flow. 

B. Refermionization 

We now complete our effective single-mode description 
by refermionizing the theory (88). For that we first present 
this Lagrangian in the canonical form by rescaUng the charge 
mode: 

&^29° = ELi©- . (89) 

The field 6 is the canonical displacement field for the total 
density, 

P = -d,e . (90) 

TT 

To preserve correct commutation relations, the conjugate mo- 
mentum lie for the field (89) should be half of that for the 
field e°. 

Us = in^o = i eLi na , (91) 

with n„ defined in Eq. (27). Changing variables in the La- 
grangian (88) according to (89) and (91) we obtain the effec- 
tive Lagrangian for the mode Q, 

-i/ifl'cos(2e -I- mtot^exta; - 2A' sinke^^x) > 
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with the redefined parameters 

v' = 4v, K'=^K, A' = AA 



A 



K'hke^iv' ■ ^^^^ 

The rescaling (93) has the following meaning. The velocity 
quadrupling simply states that we count incoming fermions 
regardless of their original flavors. The rescahng of the charge 
stiffness can be understood by noting that, with the same accu- 
racy that has allowed us to discard the flavor modes [K 1, 
see above], 



1 



(94) 



which is by definition the charge stiffness for the spinless 
Dirac fermions of velocity v' and density of states p' [cf. 
Eq. (33)]. Finally, the external fields U and /x are not rescaled 
[cf. Eqs. (39) and (37)]. 

Refermionizing the Lagrangian (92) by introducing the 
Dirac spinors ^ ~ (27ra')~^/^e'®, we formaUy obtain the 
effective Hamiltonian 

Weff [*] = Jdx {-ihv'a^d^ + AVi + U{x) 

-m} * + J I] P~kV{k)pk ^^^^ 



in which case the charge excitation is extended over several 
potential minima, and the single-mode phase soliton approach 
[Sec. V C] applies, with the system (92) of the form (54) with 
K ^ K',b 2A', 13 = 2, and wifli the period ratio defined 
in terms of the total density, m mtot- 

The excitation gaps are then obtained in a standard way: 
According to Sec. V C, when the total density mtot is integer, 
one averages the potential term in Eq. (92) over the period 
obtaining the effective coupUng 



9 = gJm,A2A'), 

then integrates over the fluctuations of the 6-field 



(100) 



g{l)=g{a'){- 



'K' 



4?7«:l, (101) 



with K' given by Eq. (94) and rj by Eq. (68), finds the corre- 
sponding charge soliton scale self-consistently as 



(102) 



and obtains the excitation gap 2Am, where 



for the fictitious spinless Dirac fermions. These are the NT 
electrons averaged (traced) over their SU(4) flavor configura- 
tions with equal weights. Naturally, the total fermion number 
density (6) in the new variables 



p{x) = vp+vp , 



(96) 



The effective gap A' in Eq. (95) is chosen in such a way that 
it corresponds to the renormalized coupUng (45) entering the 
Lagrangian (92), 



/i 



A' 



-5- 



a' ~ Zfl 



(97) 



similarly to the definition (38). In Eq. (97) the new length 
scale cutoff a' is assumed since the present approach is valid 
only at the length scales I > la beyond the correlation length 
of the neutral sector.^^ Eq. (97) together with Eqs. (45) and 
(64) estimates 



A' ^ ^ ^ D'/'At^' . 



(98) 



To summarize, the NT electron dynamics in the dilute limit 
(61) is described by the effective Hamiltonian (95) of spin- 
less Dirac electrons. These fictitious fermions of the density 
(6) carry unit charge, and interact with each other and with 
external fields in a standard way. 



C. Excitation gaps 

The (charge) excitation gaps can be now estimated from 
the effective Lagrangian (92) or from the Hamiltonian (95). 
Below we assume sufficiently long range interaction 



^ch > Ag 



(99) 



Here the charge gap 



A, mtot = 4m = ±1, ±2,.... 

(103) 



(104) 



with A' given by Eq. (98), and Jmtot = Jim is the Bessel 
function that depends on the total density (1). 

When the density mtot = p' /q' is a. simple fraction, 
the charge density in the commensurate configuration is 
(Z'Aext— periodic [as shown in Fig. 1(a) for mtot = 1/2]. The 
gaps then foUow from the standard phase soUton approach^^ 
[Sec. V C] which we will not repeat here. 



D. Classical limit 

The classical limit of the Hamiltonian (95) is obtained by dis- 
carding the 0{rj') quantum fluctuations, and coarse-graining 
beyond the length scale ^a, effectively treating electrons as 
distinguishable point particles [on the scale I > Ir] that inter- 
act with each other via the Coulomb potential (7). The classi- 
cal energy of this system 



+ Y^v{xi-Xi,) + Y,v{y3 



- Y.[A'-U{y,)] 
y,)-Y.yix,-y,)M^'^ 



Here the indices i, i' and run over electrons (with posi- 
tions x) and holes (y), correspondingly, in the minima and in 
the maxima of the potential (14). The incompressible states 
are parametrized by the pairs (ug, rih) of rational numbers of 
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electrons and holes per period Aext according to Eq. (2) [cf. 
Section 1 and Fig. 1]. 

The connection between the classical limit (105) of the ef- 
fective Hamiltonian (95), and the original charge mode La- 
grangian (92) follows from representing the coordinates of 
electrons by 



— J- wi, 



(106) 



where (j)j are the displacements from the ideal positions 



(0) _ jAext 



mtot 



(107) 



in the Wigner crystal in the absence of the external poten- 
tial. [For simplicity, we are not considering the holes in the 
maxima of U{x), which can be treated analogously.] In the 
continual limit = (p{xj) « ^{x) the charge density is 



p{x) 



(108) 



describing the change of (j> by Aext / "^tot when one extra par- 
ticle is added. The third sum in Eq. (105) gives the interaction 
energy between the electrons. Substituting the coordinates 
(106) into this sum and expanding it up to the second order 



„(o) 



iiKpi — pa (xf^^ — Xj'' j dx(j), we obtain the gradient term 
for the displacement mode energy, 

\ jdx (2e^lniv)p^(.)^^ jdx^^^, (109) 



where 



e = — — (f>{x) 

/^ext 



(110) 



is the net charge mode (89) defined in accord with the expres- 
sion (90) for the total charge density, and the charge stiffness 



k' 



1 



(Ill) 



Here the stiffness K' is defined in Eq. (94), and the argument 

TV of the Coulomb logarithm is found self-consistently as a 
number of electrons whose positions are altered in a charged 
phase-soliton excitation. The expression (109) is the classical 
limit of the second term of the Lagrangian (92) with a meaning 
of the Coulomb interaction between the quasi-classical elec- 
trons. The nonlinear term of the Lagrangian (92) gives the 
bare energy cost A' cos (rritot^ext^^) = ^' of adding a parti- 
cle above the gap (modulo interactions), with Scf) = Aext /mtot 
according to Eq. (108). This term corresponds to the constant 
term in the classical energy (105), A'. The interaction 
with the potential, U{xi), can be written in the usual 
form (13) and added to the argument of the cosine via the 
gauge transformation (32). 



E. Discussion 

The main result of this Section is the single-mode effective 
Hamiltonian (95) that has allowed us to map the problem of 
the interacting fermions with the four flavors onto that of a 
single flavor and to utilize the standard phase soliton approach 
for the single mode. 

The physical meaning of the present treatment is as follows. 
In the noninteracting case, fermions of the same flavor avoid 
each other due to the Pauli principle. The ground state wave 
function is then given by a product of the four Slater determi- 
nants, one for each flavor. However, when the repulsion be- 
tween the fermions is strong, fermions of all the flavors avoid 
each other in a similar way, and the ground state wave function 
is a Slater determinant of a four-fold size.^^ This is manifest 
in the wtot = 4m— dependence of the excitation gaps (103). 

The original SU(4) flavor symmetry of the problem be- 
comes manifest on the level of renormalization, namely in the 
particular scaling law of 4/5 of the renormalized gap (98) pro- 
duced by the renormalization group flow of the flavor sector 
on the length scales a < Z < Zfl. 

Let us compare the scaling laws in Eqs. (103) and (104) 
with those obtained in the weak-coupling limit. In the classi- 
cal limit the charge excitation gap (104) for the stand-alone 
nanotube has the same form as that obtained in Sec. VE 
above, Eq. (67): 



A = A when K ^ oo . 



(112) 



However, one notes that at finite K, the power law exponents 
in the expressions (104) and (103) are different from those 
in Eqs. (71) and (73). This is not surprising since the the- 
ory (88) is different from the original model (35). Whereas in 
the latter the exchange is important, in the former the flavor 
sector is traced over under the assumptions specified above 
in Sec. VII A; different theories indeed yield different scaling 
laws. A similar subtlety in taking the order of Umits of van- 
ishing temperature and exchange was observed in the recent 
calculation of correlation functions^^ for spin— 1/2 fermions, 
and was qualitatively explained in Ref . 44. In the present case, 
the same phenomenon is manifest in the scaling behavior of 
the excitation energy. 

We underline that the external potential (14), by adding an 
extra length scale Aext, naturally distinguishes the regimes 
Zfl > Aext and Za < Aext in which the flavor physics is im- 
portant and unimportant, correspondingly. Accordingly, the 
dependence on the parameters of the potential in these two 
regimes is also qualitatively different, as one may see by com- 
paring e.g. Eqs. (73) and (103) even in the limit K ^ oo. 



Vni. PHASE DIAGRAM 

Below we draw the phase diagram for the classical system 
(105) in the A, p plane. Fig. 6, where A is the potential am- 
plitude in Eq. (14). Here we assume that the screening length 
is larger than the length L of the system. 



Zs ~ /ch > i = iVAe 



(113) 
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(effectively meaning that Ich = oo), thus the charged phase- 
soliton optimization of an excitation (described in Sec. V C) 
does not take place. In the limit (113), by raising the chemical 
potential one adds electrons (holes) to the sequence of "quan- 
tum dots" ( "anti-dots") evenly over all the system, chang- 
ing the commensurate density mtot from one rational number 
to another Each region in Fig. 6 corresponds to a particular 
commensurate phase (rig, Uh) (as described in Section 1), with 
density (2). Borders separating regions with integer (rie, Uh) 
are comprised of fractional-density states, such as the domain 
(1/2,0) between (0,0) and (1,0) [charge state sketched in 
Fig. 1(a)]. At zero temperature, all the commensurate charge 
states (rational mtot = p' / <l') are incompressible, with a spec- 
trum being a devil's staircase. At finite temperature gaps for 
sufficiently large denominator q' will be washed out. The un- 
derlying Dirac symmetry makes the phase diagram symmetric 
with respect to /i — /x, rig ^ rih, so that only the /i > 
part is shown. 

The phase diagram is obtained by minimizing the energy 
functional 



(114) 



[Ec\ given by Eq. (105)], with respect to positions Xi and yi of 
electrons and holes in the following way. For large size N, we 
neglect finite size effects and utiUze translational invariance. 
In this case, in a commensurate state, the optimal positions of 
electrons and holes relative to each potential minimum (maxi- 
mum) are the same for each period. The result is the family of 
system's charging energy values £'„^,„^(A) per period Aoxt, 
such as 



E 



10 



A' -A 



e 

Acxt 



N/2 

^ n 

n—l 



Eu 



2Eif) — 



2e^ 



N/2 
n=l 



(115) 



(116) 



The first term in (1 18) is the energy of Ue electrons placed into 
each minimum and Uh holes into each maximum of U (x). It 
corresponds to the first two terms of (105). The second term 
in (118) is the interaction energy of the electrons and holes 
located in the extrema of U (x). Here 



Co 



A. 



21n(i/Aext) 



(119) 



is the NT capacitance per period Acxt- Finally, in (1 18) is 
the interaction energy of n electrons (or n holes) minimized 
with respect to their positions inside the corresponding poten- 
tial well of the periodic potential (14). The Dirac symmetry 



yields En 



— En 



Minimization of Eq. (114) using 



the approximation (118) yields the phase diagram sketched in 
Fig. 7. In the limit (113), excitation gaps corresponding to the 
incompressible states with the total density n = rie — Uh = 
(rie + 1) — {uh + 1) = ... = + s) — {uh + s) oscillate but 
do not vanish. Their minimum value 



e 

C'o 



(120) 



is determined by the NT charging energy. The regions 
(rie, Tih) and (rig + 1, rih + 1) of the phase diagram are sepa- 
rated by the vertical hnes of fixed A, with its value implicitly 
determined from 

^ = A' + i {K. + l(A) - VnM) + Vn, + l{A) ~ VnAA)} . 

(121) 

Away from the values (121) the gap increases (x A. 

The mean field approximation (118) produces fairly accu- 
rate borders between the domains of the phase diagram in 
Fig. 6. Below we consider the quantum-dot addition energies 
Vn^ for the small and large rig. 



A. 



N/2 

E 



E20 = niin < 2 A' — 2 A cos k^xtSx + 



2Sa 



26x/K 



— 26x/Xc 



(117) 



[here zLSx is the electron coordinate relative to the minimum 
of U{x)], and so forth. While minimizing the energy En^.nh 
with respect to the positions of electrons and holes within each 
potential period, unphysical configurations of electrons on top 
of the holes are excluded by demanding that their minimum 
separation be e^/A'. The latter condition takes into account 
exciton binding energy in the Dirac system. 

The functional En^.n^ {A) can be further approximated by 
minimizing the interaction energy between the charges only in 
the same potential minimum (or maximum), treating the rest 
of the system in a mean-field way: 

(rig - nh)^e'^ 



En,,n^ w {ne + ni,){A' ^A) + - 



2Cn 



(118) 




FIG. 6: (Color online) Phase diagram for the nanotube, A' = 
67reVAext, ln(L/Aext) = 2.5, e^Aext = 1.44 meV for Acxt = 
l/im. Shaded regions are the result of the numerical minimization 
of Eq. (114) (see text), black lines are approximations obtained from 
Eq. (118) with V„ given by Eq. (123). 
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FIG. 7: Phase diagram in thie (A, ji) plane according to the model 
(118). The minimum width of the gap for the incompressible state 
with any integer density rie — nn is given by the NT charging energy 
(120). In the limit Tie, 3> 1 all the regions of the phase diagram 
are identical, with their period in A given by Eq. (128). 



A. Small density ne,nh ~ 1 

First consider the situation when the Coulomb interaction 
over a period is small, e^/Aext ^ ^> such that rig electrons 
are grouped close to the minima of the potential that can be 
approximated by a piecewise-quadratic polynomial: 

U{x) « -A + min„<iv^^ {x-in+ i)Aext) • (122) 

Minimizing the Coulomb energy of rie charges in a single 
minimum of the potential (122), one obtains the values 



Vo = Vi = 0, 



V3 = 5y^V2iA) , ... . 



(123) 



The power law correction ^ A^^"^ due to the above expres- 
sions is observed in Fig. 6 for A ~ A' as a deviation from 
the straight lines that separate different regions of the phase 
diagram. 

When, on the other hand, the periodic potential is a small 
perturbation that shifts electrons from their equidistant posi- 
tion in a Wigner crystal (A <C e^/Aext)> the perturbation the- 
ory in A yields analytic behavior of the gap widths on A. 



B. Lat^e density ne,nh » 1 

At large rig, :» 1 the model (118) can be simplified by us- 
ing the continuous Thomas-Fermi description for the density 
of the classical electrons (holes) inside each potential min- 
imum (maximum). To find Vn^, we approximate the local 
charge density as 



p{x) w ptf = -T — cos fcexta; , 

Aext 



(124) 



that mimics the external potential profile (14). The density 
(124) is normalized to Ue = J^^''^J^'^ pTF{x)dx. The inter- 
action energy Vn^ or Vn^ inside each "quantum dot" can be 
estimated as 



Aext/4 \X-X'\ 



2 7 

2cr 



where the "dot capacitance" 



27r In- 



(125) 



(126) 



is independent of the potential amplitude A. The log singular- 
ity in the integral in Eq. (125) is cut off on the scale a' ^ Zfl 
of the order of the flavor soliton size, below which the quasi- 
classical description breaks down (see Section VII). The A- 
dependence of Ci would appear as a correction to y„_. with 
^-dependent integration hmits in Eq. (125). The case of 
the holes in the potential maxima ("anti-dots") is analogous. 
Eq. (121) then yields the potential amplitude values that sep- 
arate the configurations (ug, rih) and (rig + 1, ri/i + 1): 



i(ne>"fc)->("e + l>")i + l) 



A' + — K + n^ + l). (127) 



In this limit aU the honeycomb regions in Fig. 7 are identical 
with their size in A being 



5Ar^ 



e 



(128) 



Qualitatively, the regular honeycomb structure of Fig. 7 ap- 
pears already for moderate amplitudes A (Fig. 6). The bor- 
ders between the regions of the diagram in the limit A > 
e^/Acxt, A', are approximately linear, dominated by the lin- 



ear dependence of En^^nn 
ofEq. (118). 



on A that stems from the first term 



C. Crossover with tlie case kh < L 

We will now show that the asymptotic values of the positions 
(127) of the gap minima derived above for large and Uh 
coincide with those obtained in Section VII utilizing the phase 
soUton method in an infinite system. For this crossover we 
assume Ich ~ ~ Aext- 

Consider the positions of the minima for excitation gaps 
(103) (as a function of the potential amplitude) 



i' = 4i = 



4A 



TT 

2 



(129) 



where we took the value of the charge stiffness (33) at the 

momentum scale A^^^ and substituted K^K — 1 at 

1. Utilizing the asymptotic behavior of the n— th zero of the 

Bessel function Jm^„^ , 



mtot 



const 



Trmtot 



7rn, n = min {rie, n/i} , 
(130) 

we indeed obtain the values A ~ {n^ + n^)e^/2Ci corre- 
sponding to gap minima in accord with Eq. (127). 
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K. EXPERIMENTAL MEANS TO PROBE 
INCOMPRESSIBLE STATES 

A. Conductance measurements 

Excitation gaps corresponding to the incompressible states 
can be detected in transport measurements. When, by means 
of varying the gate voltage, the Fermi level is tuned to be in 
the gap, the zero-bias conductance a across the nanotube van- 
ishes. At finite temperature T, u will have the activated form 

<T ~ e-'^»'=*/'^ , (131) 

with the activation energy Tact = equal to half of the 
excitation gap corresponding to a particular incompressible 
state. Physically, Eq. (131) corresponds to the transport by 
means of thermally-excited charge solitons that can carry elec- 
tric current through the nanotube. (A similar situation has 
been recently considered in the context of transport in gran- 
ular arrays. ^^) The values and positions of the gaps can be 
also revealed by measuring the differential conductance at fi- 
nite bias with varying gate voltage. In this case, peaks in the 
differential conductance as a function of the applied bias will 
mark the positions of the excitation gaps. 

B. Adiabatic charge pumping 

A more challenging possibility is to realize the Thouless 
pump^^ in a nanotube. Such a setup requires an adiabatically 
moving periodic potential that could be created e.g. by cou- 
pling to a surface acoustic wave (SAW),^'^ or by sequentially 
modulated gate voltages on the array of underlying gates. 
When the chemical potential is inside the rn-th minigap A„i, 
the adiabatically moving periodic potential with frequency / 
induces quantized currenf^ 

3 = mtotef , mtot = 4m . (132) 

In such a setting, novel fractional-m incompressible states 
considered in this work will manifest themselves in the adi- 
abatic current (132) quantized in the corresponding/racf/o«.s 
of 4e/.^^ This result can be understood by invoking the topo- 
logical invariant property of the Thouless current.'^^''*^ The ex- 
pression (132) is trivial in a semiclassical limit, with a mean- 
ing of transporting on average mtot electrons per cycle in a 
conveyer-belt fashion. By staying inside the gap and adiabati- 
cally changing the parameters of the system, the current (132) 
remains invariant and hence it is valid in the fully quantum- 
mechanical case. 

Operation of the charge pump requires adiabaticity 

ksT, hf<^Am. (133) 

Since the energy scale for the minigaps is set either by the 
gap 2Ao at half-filhng, or by the strength of the Coulomb in- 
teractions between electrons separated by ~ Aext> the typi- 
cal minigap values A„j are in the meV range, and the adi- 
abaticity condition (133) is realistic. The feasibility of the 



Thouless pump in the NT-SAW setup is further corroborated 
by recent pumping experiments involving SAWs. In particu- 
lar, in the pumping of electrons between the two 2D electron 
gases through a pinched point contact"*' the achieved quality 
of current quantization is close to metrological.^'''^' Recently 
the SAW-assisted pumping has been demonsttated through the 
laterally defined quantum dot,^^ as well as through the semi- 
conducting nanotube whose working length L matched the 
SAWperiod, L = Xcxt-^^'^'* 

We contrast the non-dissipative current (132) on the plateau 
with the dissipative non-quantized current away from com- 
mensuration considered in Ref. 55 for arbitrary quantum wire. 
Such a non-quantized current is pumped when no gap opens 
(for incommensurate densities or for commensurate densities 
with the interaction below criticality) and is characterized by 
the interaction-dependent critical exponents. 

A practical realization of the proposed pumping setup can 
become a first implementation of the Thouless transport. Be- 
sides being instrumental in studying electron interactions in 
a nanotube (by detecting and measuring fractional— m mini- 
gaps that arise solely due to interactions), it could realize the 
"conveyer belt" for electrons with a possibility of pumping 
current quantized in fractions of the unit charge per cycle. 
In other words, such a setup would make the first example 
of the charge pump operating at the. fraction of the base fre- 
quency. This setup could allow one to study in detail the elec- 
tron correlations and interactions in nanotubes, with a variety 
of controllable parameters at hand, such as the shape and fre- 
quency of the external potential, nanotube gap Ag (that can be 
modified externally), and the gate voltage. By exploring the 
phase diagram one can study effects of Wigner crystallization, 
quantum commensurate-incommensurate transitions, and the 
Tomonaga-Luttinger correlations. 

Finally, we note that the described setup can be utilized to 
adiabatically transport low-energy strongly correlated SU(4) 
flavor states [e.g. those for m = 1/2 obeying the effec- 
tive Gross-Neveu Lagrangian of the form (75) described in 
Sec. VI] over a macroscopic distance, since the coupling to the 
adiabatically moving external potential is SU(4) invariant and 
thus it does not desttoy spin or flavor correlations. This may 
become useful in the context of solid state implementations of 
quantum information processing. Furthermore, this setup can 
provide a possibility of realizing a quantized spin-polarized 
pump by subjecting the system to magnetic field tuned in such 
a way that only one spin population of electrons is commensu- 
rate with the potential and participates in current quantization. 



X. CONCLUSIONS 

In the present work we have shown that coupUng of the 
interacting electrons in a nanotube to an external periodic po- 
tential is a rich setup to study one-dimensional electronic cor- 
relations, including the crossover between the Luttinger liquid 
and the Wigner crystal. We demonstrated that the external po- 
tential locks the system into incompressible states. The corre- 
sponding excitation gaps (estimated to be in meV range) are 
found by adequately treating the curvature of the electronic 
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dispersion in the bosonized language, and by further gener- 
alizing the phase soliton method onto the case of multiple 
modes. 

In the regime when the gaps open due to the Bragg diffrac- 
tion in a multi-flavor Luttinger Uquid, we identified and inves- 
tigated the novel incompressible fractional-density state with 
m = 1/2 electrons of each flavor per period of the potential. 
The phase soliton action derived for this case describes the 
charge excitation, and the SU(4)-flavor excitations governed 
by the 0(6) Gross-Neveu model. 

In the opposite Umit of the pinned Wigner crystal we de- 
rived the effective single-mode Hamiltonian and found that 
the phase diagram in the classical limit has a stucture of the 
devil's staircase. 

The interaction-induced incompressible states can be de- 
tected in the Thouless pump setup that can allow one to study 
electron correlations and the transition to the Wigner crystal, 
as weU as to realize the quantized charge pump that operates at 
a fraction of the base frequency by virtue of electron-electron 
interactions. 
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APPENDIX A: CONDITIONS FOR THE TIGHT-BINDING 
LIMIT 

Below we derive the condition (26) for the tight binding limit 

(25) in the single particle picture. We consider the two cases 
depending on the relation between the potential amplitude A 
and the Dirac mass term Aq. 

The case A <C Aq corresponds to the "nonrelativistic limit" 
of the Dirac equation, where all the relevant energies are much 
smaller than the gap Aq. We define the "Dirac" electron mass 
M via Aq = Mv'^, and turn to the effective Schrodinger de- 
scription with the Hamiltonian 



ch 



2M 



+ U{x) , 



(Al) 



in full analogy with the non-relativistic limit in QED. The tun- 
neUng amplitude between the adjacent minima of U [x) is pro- 
portional to e~^^^<^o/^^ where the classical action under the 
barrier 



'S'U<Ao 



(A2) 



Therefore both the minibands and tunneling are suppressed in 
the limit A < Aq if 



Ao 
eo 



A Ao 
< — < — 
eo eo 



The inherently Dirac regime occurs when A > Aq. In this 
case electron can tunnel between the minima of the potential 
(14) sequentially through the hole part of the spectrum. The 
corresponding tunneUng ampUtude ~ e-25'U>Ao/?i^ where 
the classical action for the particle with energy E under bar- 
rier between the electron and hole parts is 

S\a>a, = \jdx ^J\[E-U{x)f-Al\. (A4) 

Tunneling from the minimum of U{x) at E = Aq — A to the 
hole part of the spectrum yields (6 = Aq/A) 




(A5) 



Tunneling from the energy level \E\ ^ A — Aq far 
from the potential bottom yields the action <S'||£;|<a-Ao — 
ttAq/UvS, where 5 = \dxU{x)\ ~ fcext^- The requirement 
S\e=Ao-a > h or S'j|£;|^4_Ao > ^ yields that the single 
particle bandwidth in the Dirac regime is exponentially sup- 
pressed if 



^ < ^ < 
eo eo V eo 



(A6) 



where, correspondingly, the exponent c = 3 for tunneling 
from the potential minimum and c = 2 for tunneling from 
energy level far from the potential bottom. 

Summarizing, Eqs. (A3) and (A6) yield the condition (26) 
for the tight-binding Unlit (25). 
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PHASE SOLITON ACTION 



Below we consider the weak coupling limit for the fractional 
density m = i, characterized by the chemical potential /i = j 
in accord with Eq. (40). 

Our course of action has been outlined in Section V above. 
Technically, it will be more convenient to first work in the ba- 
sis of the original fields <da and later utilize the transformation 
(29) to obtain the Lagrangian £1/2 ) ^"l in the charge-flavor 
basis. Since the average of the potential energy (37) over the 
two successive potential periods is zero, we decompose the 
fields Qa (x) in a series in the coupling go keeping the zeroth 
and the first orders 



= + e^i) , ei^) = o{go) . 

Consider the Hamiltonian 

K -1 



(Bl) 



(A3) +^go cos(2ec + ^h^tx - 2Asm fcexta;)| • (B2) 
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The Euler-Lagrange equations (57^^^/^^ /^Oq = are 



K - 1 

a. XX ' ~A ^xx 



e(i) 



5o . 



s\-B.{2QoL+\K:^tX-2A sin fccxt^;), 

(B3) 

where 5* = ^ 9q, and Sx = dxS, etc. Integrating Eq. (B3) 
we obtain 



S'(i) - 



1 - K 



4 +A:e.t 
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J„i(2yl) 

2m 



(B4) 
(B5) 



cos 



(2e„ + (i-m)fcexta;).(B6) 



Substituting Eqs. (B4) and (B6) into the Hamiltonian (B2), 
after somewhat lengthy but straightforward algebra the slow 
mode potential follows: 



^1/2 = 



IBtt 



(4- K)i;i/2^cos4ea 



+ '*ui/2 E cos(26q — 28q') 

— K,Vi/2 E C0s(2G)a + 29q' 



(B7) 



Here 



ffo 



.90 

K- 1 

if 



Ao 
eoa 



and the couplings Vi/2 and Ui/2 are defined as 

Jm(^)Jl-m(^) 



Vl/2{z) = E 



(2to- 1)2 



m— — 00 



(B8) 
(B9) 
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FIG. 8: (Color online) The functions 111/2 and Ui/2 defined in 
Eqs. (BIO) and (Bll). We find that Ui/2(z) < u^/2{z) holds for 
all z. Zeroes of V1/2: 2 = 0, ±2.33, ±3.80, ±5.47, .... 

with z being a shorthand for 2 A. The functions wi/2(z) and 
Mi/2(z) are plotted in Fig. 8. 

In the commensurate state the minimum value 



min\/i/2[e«]--^{4K/2(2i)| 



"i/2(2i)} 

(B12) 

of the potential (B7) corresponds to {Qa] being a permutation 

of a set 01 (/)2 02}, with 0i,2 = ±7r/4 for vi/2{2A) > 
and 01 = 0, 02 = 7r/2 for t;i/2 {2 A) < 0. 

We now discuss the obtained commensurate classical state. 
The ground state degeneracy in the noninteracting case (k = 
0) is equal to 2'^. This follows from the potential (B7) in which 
only the first term is nonzero. In the presence of interactions 
(k > 0), the other terms in (B7) reduce this degeneracy from 
16 to six. This result could have been foreseen without any 
calculation since the remaining degeneracy is a number of 
configurations in which any two different fields, Qa and Qa' 
with a 7^ a', are placed in the same minimum of the external 
potential. In other words, in the bosonized treatment the ef- 
fect of fermionic exchange manifests itself in a stronger repul- 
sion between the solitons of the same flavor This is in agree- 
ment with the numerical minimization performed in Sec. V B, 
Fig. 4. The symmetry of the obtained classical ground state is 
schematically illustrated in Fig. 1(e). 

The change of variables (29) in the potential (B7) yields 
Eq. (77). 



* Electronic address: dima@alum.mit.edu 
' S. lijima, Nature 354, 56-58 (1991). 

^ R. Saito, G. Dresselhaus and M. S. Dresselhaus, Physical Proper- 
ties of Carbon Nanotubes, Imperial College Press, London, 1998. 

^ M.S. Dresselhaus, G. Dresselhaus, Ph. Avouris, Carbon Nan- 
otubes: Synthesis, Structure, Properties and Applications, 
Springer Verlag, New York, 2001. 

" PL. McEuen, Nature 393, 15 (1998). 

^ C. Dekker, Physics Today, May 1999, p. 22. 

^ Y.A. Krotov, D.-H. Lee, S.G. Louie, Phys. Rev. Lett. 78, 4245 
(1997). 

^ R. Egger, A. O. Gogolin, Phys. Rev. Lett. 79, 5082 (1997). 
C. Kane, L. Balents, M. P A. Fisher, Phys. Rev. Lett. 79, 5086 



(1997). 

' H. Yoshioka and A. A. Odintsov, Phys. Rev. Lett. 82, 374 (1999). 
L.S. Levitov and A.M. Tsvelik, Phys. Rev. Lett. 90, 016401 
(2003). 

" M. Stone (Ed.), Bosonization, World Scientific, Singapore (1994) 
M. Bockrath, D. H. Cobden, P L. McEuen, N. G. Chopra, A. Zettl, 
A. Thess, R. E. Smalley, Science 275, 1922 (1997). 
S. Tans, M.H. Devoret, R. Groeneveld and C. Dekker, Nature 394 
761 (1998). 

M. Bocki-ath, D. H. Cobden, Jia Lu, A. G. Rinzler, R. E. Smalley, 

L. Balents, P L. McEuen, Nature 397, 598 (1999). 

Z. Yao, H. W. Ch. Postma, L. Balents, C. Dekker, Nature 402, 273 

(1999). 



21 



" J. Nygard, D.H. Cobden, M. Bockrath, P.L. McEuen, RE. Linde- 

lof, Applied Physics A 69, 297 (1999). 

H. Ishii, H. Kataura, H. Shiozawa, H. Yoshioka, H. Otsubo, 
Y. Takayama, T. Miyahara, S. Suzuki, Y. Achiba, M. Nakatake, 
T. Narimura, M. Higashiguchi, K. Shimada, H. Namatame and 
M. Taniguchi, Nature (London) 426, 540 (2003). 
P. Jarillo-Herrero, J. Kong, H.S.J, van der Zant, C. Dekker, 
L.P. Kouwenhoven, and S. De Franceschi, Nature (London) 434, 
484 (2005). 

M.-S. Choi, R. Lopez, and R. Aguado, Phys. Rev. Lett. 95, 067204 
(2005). 

^ V.I. Talyanskii, D.S. Novikov, B.D. Simons, L.S. Levitov, Phys. 
Rev. Lett. 87, 276802 (2001). 

^' D.S. Novikov, Phys. Rev. Lett. 95, 066401 (2005). 

22 D. J. Thouless, Phys. Rev. B 27, 6083 (1983). 

2^ An interesting quasi- Id analogy of QHE has been proposed re- 
cently by C.L. Kane, R. Mukhopadhyay, and T.C. Lubensky, Phys. 
Rev. Lett. 88, 036401 (2002). 

^ J. Frenkel and T.A. Kontorova, Zh. Eksp. Teor. Fiz. 8, 1340 
(1938). 

2^ F.C. Frank and J.H. van der Merwe, Proc. Roy. Soc. London, Ser. 

A, 198, 205 (1949). 
2* I.E. Dzyaloshinsky, Zh.E.T.R 47, 1420 (1964). 
" V.L. Pokrovsky, A.L. Talapov, Sov. Phys. JETP 75, 1151 (1978). 
^ P. Bak, Rep. Prog. Phys. 45, 587 (1982). 

2' Yu.V. Nazarov, D.V. Averin, Phys. Rev. Lett. 81, 653 (1998); 
R. Klesse, A. Stern, Phys. Rev. B 62, 16912 (2000); M. Pustil- 
nik, E.G. Mishchenko, L.I. Glazman, and A.V. Andreev, Phys. 
Rev Lett. 91, 126805 (2003). 
S. Coleman, Phys. Rev. D 11, 2088 (1975). 

" RD.M. Haldane, J. Phys. A 15, 507 (1982). 

'2 To be precise, this statement is true in the forward scattering ap- 
proximation that holds down to undetectably small energy scale, 
below which further interaction-induced gaps may open even in 
the absence of the bare gap, Ao = 0.^' 
C.L. Kane and E.J. Mele, Phys. Rev. Lett. 78, 1932 (1997). 

C. Zhou, J. Kong, and H. Dai, Phys. Rev. Lett. 84, 5604 (2000); 
M. Ouyang, J.L. Huang, C.L. Cheung, CM. Lieber, Science 292, 
5517 (2001). 

" H. Ajiki and T. Ando, J. Phys. Soc. Jpn. 65, 505 (1996). 

J. -O. Lee, J. R. Kim, J. J. Kim, J. Kim, N. Kim, J. W. Park, and 
K. H. Yoo, Sol. Stat. Comm. 115, 467 (2000). 

" D.S. Novikov and L.S. Levitov, cond-mat/0204499 (unpub- 
lished); cond-mat / 0512059 (Phys. Rev. Lett., in press). 

D. S. Novikov, unpubhshed 



K. Wilson, Phys. Rev. 179, 1499 (1969); C.G. Callan, Phys. Rev. 
D 2, 1541 (1970); K. Symanzik, Commun. Math. Phys. 18, 227 
(1970); S. Coleman, Ann. Phys. (N.Y.) 101, 239 (1976). 
Similar expression for renormalized gap A has been obtained 
earUer in Ref. 10. Whereas we agree on scaling exponents of 
D and Ao, the if -dependence is different. The procedure of 
Ref. 10 yields the prefactor K^^^ instead of the more natural re- 
sult /('I' 2+C'{»)) jj^jjj matches the classical limit. We believe this 
discrepancy is a consequence of oversimplifying the renormaliza- 
tion group treatment by introducing a single soUton scale under 
the RG in Ref. 10. 

D.J. Gross and A. Neveu, Phys. Rev. D 10, 3235 (1974). 

A.B. Zamolodchikov and Al.B. Zamolodchikov, Ann. Phys. 

(N.Y.) 120, 253 (1979). 

V.V. Cheianov and M.B. Zvonarev, Phys. Rev. Lett. 92, 176401 
(2004); J. Phys. A 37, 2261 (2004). 

G.A. Fiete and L. Balents, Phys. Rev. Lett. 93, 226401 (2004). 
Here we do not distinguish between ln{ls/a) and ln(/s/a') in the 
Coulomb formfactor V{q), with K :$> 1 implying large screening 
length Is ~ Icb ^ a.a'. This corresponds to treating the tube 
capacitance with logarithmic accuracy. 

A similar situation occurs in the Bethe-ansatz solution of the Hub- 
bard model for spin— 1/2 fermions as demonstrated by M. Ogata 
and H. Shiba, Phys. Rev. B 41, 2326 (1990). 
'^^ A. Altland, L. I. Glazman, and A. Kamenev, Phys. Rev. Lett. 92, 
026801 (2004). 

Q. Niu and D.J. Thouless, J. Phys. A 17, 2453 (1984). 
J.M. Shilton, V.I. Talyanskii, M. Pepper, D.A. Ritchie, 
J.E.R Frost, C.J.B. Ford, C.G. Smith, G.A.C. Jones, J. Phys. 
Condens. Matter 8, L531 (1996); V.I. Talyanskii, J.M. Shilton, 
M. Pepper, C.G. Smith, C.J.B. Ford, E.H. Linfield, D.A. Ritchie, 
G.A.C. Jones, Phys. Rev B 56, 15180 (1997). 

^° Q. Niu, Phys. Rev. Lett. 64, 1812 (1990). 

J. Cunningham, V.I. Talyanskii, J.M. Shilton, M. Pepper, A. Kris- 
tensen, RE. Lindelof, J. Low Temp. Phys. 118, 555 (2000). 

" J. Ebbecke, N.E. Fletcher, T.J.B.M. Janssen, F.J. Ahlers, M. Pep- 
per, H.E. Beere, and D.A. Ritchie, Appl. Phys. Lett. 84, 4319 
(2004). 

J. Ebbecke, C.J. Strobl, and A. Wixforth, Phys. Rev. B 70, 233401 

(2004) . 

^ P. Leek, M. Buitelaar, V. Talyanskii, C. Smith, D. Anderson, 
G. Jones, J. Wei, and D. Cobden, Phys. Rev. Lett. 95, 256802 

(2005) . 

R. Citro, N. Andrei, and Q. Niu, Phys. Rev. B 68, 165312 (2003). 



